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Abstract 

A self-consistent Ansatz for a new sphaleron of SU(3) Yang-Mills-Higgs theory is 



■<^j- , presented. With a single triplet of Weyl fermions added, there exists, most likely, one 

pair of fermion zero modes, which is known to give rise to the non-Abelian (Bardeen) 
anomaly as a Berry phase. The corresponding SU(3) gauge field configuration could 



take part in the nonperturbative dynamics of Quantum Chromodynamics. 
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1 Introduction 

The Hamiltonian formalism makes clear that chiral anomalies [1-7] are di- 
rectly related to particle production [8]. The particle production, in turn, 
traces back to a point in configuration space (here, the mathematical space 
of three-dimensional bosonic field configurations with finite energy) for which 
the massless Dirac Hamiltonian has zero modes. This point in configuration 
space may correspond to a so-called sphaleron (a static, but unstable, classical 
solution), provided appropriate Higgs fields are added to the theory in order 
to counterbalance the gauge field repulsion. 

The connection of the triangle (Adler-Bell-Jackiw) anomaly [1,2] and the 
sphaleron S [9-11] is well-known. An example of an anomalous but consistent 
theory is the chiral SU(2) x £7(1) gauge theory of the electroweak Standard 
Model, where the anomaly appears in the divergence of a non-gauged vector 
current (giving rise to B + L nonconservation [5], with B and L the baryon and 
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Fig. 1. Configuration space: vacua and sphalerons. 

lepton number). The electroweak sphaleron S [11,12] is believed to contribute 
to fermion-number-violating processes in the early universe. 

The point S in configuration space locates the energy barrier between vacua 
V and V, where V is the gauge transform of V with a gauge function corre- 
sponding to the generator of the homotopy group 7i3[SU(2)}. The sphaleron S 
has one fermion zero mode for each isodoublet of left-handed Weyl fermions 
[13-17]. A one-dimensional slice of configuration space is sketched in the left 
part of Fig. 1. 

The nonperturbative Sp(n) (Witten) anomaly [6] is also related to a 
sphaleron, namely the solution S* constructed in Ref. [18]. The electroweak 
sphaleron S* sets the energy barrier for the Witten anomaly (which cancels 
out in the Standard Model) and may play a role for the asymptotics of per- 
turbation theory [19] and in multiparticle dynamics [20]. 

An example of an anomalous and inconsistent theory is the chiral Sp(2) = 
SU(2) gauge theory with a single doublet of left-handed Weyl fermions. In this 
case, there is a noncontractible loop (NCL) of gauge field vacua (corresponding 
to the nontrivial element of 7r 4 [S77(2)]) at the edge of a two-dimensional disc 
in configuration space. For one "point" of the disc (S* in the Yang-Mills-Higgs 
theory), there are crossing fermionic levels. The sphaleron S* has indeed two 
fermion zero modes [20,21]. The Witten anomaly results from the Berry phase 
factor —1 for a single closed path around these degenerate fermionic levels [22]. 
This phase factor gives a Mobius bundle structure over the edge of the disc 
(the NCL of gauge field vacua), which prevents the continuous implementation 
of Gauss' law. In other words, the nontrivial Berry phase factor presents an 
insurmountable obstruction to the implementation of Gauss' law and physical 
states are no longer gauge invariant [6,8]. A two-dimensional slice of SU(2) 
Yang-Mills-Higgs configuration space is sketched in the middle part of Fig. 1. 

There remains the perturbative non-Abelian (Bardeen) anomaly [3,4,7]. An 
example of an anomalous and inconsistent theory is the chiral SU(3) gauge 
theory with a single triplet of left-handed Weyl fermions. Now, there is a non- 
contractible sphere (NCS) of gauge field vacua (corresponding to the generator 
of 7r 5 [S77(3)]) at the border of a three-ball in configuration space. By the fam- 
ily index theorem, there are crossing energy levels of the Dirac Hamiltonian 
at one "point" of the three-ball [8]. The non-Abelian anomaly is the Berry 
phase [22] from these degenerate fermionic levels. 

Again, a nonzero Berry phase (mod 27r) presents an obstruction to the im- 
plementation of Gauss' law. The Berry phase of a loop on the NCS considered 
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Fig. 2. SU (3) Bardeen anomaly as Berry phase. 

does not vanish and is, in fact, equal to one half of the solid angle subtended 
at the "point" with crossing energy levels by the loop; see Refs. [8,22] for 
further details and Fig. 2 for a sketch. It has been conjectured [23] that this 
special "point" in configuration space corresponds to a new type of sphaleron 
in SU(3) Yang-Mills-Higgs theory (denoted by S in the right part of Fig. 1) 
and the goal of the present article is to establish the basic properties of this 
sphaleron. 

The outline of our article is as follows. Section 2 recalls the action of chiral 
SU (3) gauge theory, mainly in order to establish notation. Section 3 discusses 
the topology responsible for the new sphaleron, together with approximate 
bosonic field configurations. Section 4 presents the generalized Ansatz for the 
bosonic fields, which can be shown to solve the field equations consistently. 
Section 5 gives the corresponding Ansatz for the fermionic field, together with 
analytic and numerical results which suggest the existence of two fermion 
zero modes. Section 6 presents some concluding remarks. Three appendices 
give the Ansatz energy density for the bosonic fields, an argument in favor of 
a nontrivial regular solution of the reduced field equations, and the fermion 
zero modes for a related gauge field background. 



2 Chiral SU(3) gauge theory 

In this article, we consider a relatively simple theory, namely 577(3) Yang- 
Mills-Higgs theory with a single triplet of complex scalar fields and a single 
triplet of massless left-handed Weyl fermions, both in the complex 3 represen- 
tation of S77(3). Our results can be readily extended to other gauge groups G 
with 7r^[G] = Z and appropriate representations. 1 

The chiral SU (3) gauge theory has the following classical action: 

S = J^d 4 x ^tv F^F^ + (D^y (D^) - X ($t$ - v 2 y + i^/J^} , 

1 Non-Abelian anomalies also appear in certain (2n)-dimensional chiral gauge the- 
ories with finite 7T2 n +i[G]. The requirement vr2 n +i[G] = Z is only a sufficient condi- 
tion for having non-Abelian anomalies in 2n spacetime dimensions; cf. Section 5 of 
Ref. [7]. 
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(2.1) 



where = d^A u — d u A^ + g[A^, A u ] is the SU(3) Yang-Mills field strength 
tensor and = (<9 M + g^) the covariant derivative for the 3 representation. 
The SU (3) Yang-Mills gauge field is defined as 

A,(x) = A%{x) V(2i) , (2.2) 
in terms of the eight Gell-Mann matrices 
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The field $(x) is a triplet of complex scalar fields, which acquires a vacuum 
expectation value 77 due to the Higgs potential term in the action (2.1). The 
field i/)(x) is a triplet of two-component Weyl spinors and ip(x) its complex 
conjugate. Also, we define (<x M ) = (1, —a) in terms of the Pauli matrices a 
and use the Minkowski space metric with g^ v {x) = diag(+l, — 1, — 1, —1) and 
natural units with H = c = 1. 

For static three-dimensional field configurations, we generally employ the 
spherical polar coordinates (r, 9, <p) but sometimes also the cylindrical coordi- 
nates (p, 0, z), with p = r sin 9 and z = r cos 0. 



3 Topology and approximate sphaleron 

The basic motivation for our search of a new sphaleron (S) in SU (3) Yang- 
Mills-Higgs theory has been given in the penultimate paragraph of Section 1 
and was summarized by Fig. 2. We start with the bosonic fields of the theory 
(2.1) and introduce a topologically nontrivial map S 3 A S 2 — > 577(3), with 
parameters (ip,/i,a) and coordinates (9,<f>) on the "sphere at infinity." Here, 
A denotes the so-called smash product which topologically transforms the 
Cartesian product S p x S m to S p+m by contracting, for fixed points Xq G S p 
and yo G S m , the set {x } x S m U S p x {y } to a single point. Topologically, 
one has therefore S p A S m ~ S p+m ; cf. Ref. [21]. 
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The procedure to be followed is similar to the one used for the sphalerons 
S [10,11] and S* [18], which employed maps S 1 A S 2 -> SU(2) and S 2 A 
S* 2 — > SU(2), respectively. The reason behind this particular choice of smash 
products is that the search of new solutions is best performed in a fixed gauge, 
which we take to be the radial gauge (hence, 9 and <j) appear in the maps, but 
not r). The three-ball of Fig. 2 or the right part of Fig. 1 becomes in the 
radial gauge a three-sphere (parametrized by vp, [/,, and a) with S at the top 
and V at the bottom, where the "height" corresponds to the energy of the 
configuration. See, e.g., Refs. [10,21] for further details. 

The following "strikingly simple" map [24] is a generator of n 5 [SU(3)] = Z: 



U:S 5 ^SU(3), U(z 1 ,z 2 ,z 3 ) 



( 



Z\ Z 2 + Z% 



Z\ Z2 — Z3 Z\ Z3 + Z 2 



\ 



Z\ Z3 — z 2 Z2 Z3 + Z\ 



Z2 Z% — Z\ 



(3.1) 



with Zi,z 2 ,z 3 G C and \zi\ 2 + \ z 2 \ 2 + | ^3 1 2 = 1. An appropriate parametrization 
for our purpose is 



Z2 
\Z 3 J 



/ 1 — cos 2 I ( 1 — cos tp) + i cos I sin ip cos \i ^ 

e 1<?i sin I cos | (1 — cos^) 
^ cos I sin ip sin ji (sin a + i cos a) ) 



(3.2) 



with angles ip,fi,6 G [0, ir] and a, G [0, 27r]. At first sight, Eq. (3.2) consti- 
tutes a map S 3 x 5 2 -> S 5 , with S 3 parametrized by spherical polar coordi- 
nates (ip,fj,,a) and S 2 by (9,<f>). But the expression on the right-hand side of 
Eq. (3.2) corresponds to the single point (z 1 ,z 2 ,z 3 ) = (1,0,0) if either ip — 
or 6 = 7r and the map (3.2) is effectively S 3 A S 2 — > S* 5 . The fixed points 
x G S* 3 and y G S* 2 needed for the smash product are thus given by ip = 
and 9 = Ti, respectively. 

For ip = 0, the matrix U(zi, z 2 , z 3 ) is independent of the parameters ji and 
a and the coordinates 9 and 0, 



V = U(rl>,n,a,0,<l>) 



V>=o 



1 
0-1 
0+10 



\ 



(3.3) 



For tp = ti, the matrix U (zi, z 2 , £3) is independent of \i and a but does depend 
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on 9 and 0, 

W(9,<j>) = U(il>,»,a,e,<J>) 
( cos 2 6 



1p=7T 



- cos 9 sin 9 sin e -i * N 
-cosflsinfle^ sin 2 #e 2i< ^ cos# • (3.4) 

^ -sin0e _i * -cos# , 

These two matrices have the following rotation and reflection symmetries: 

fyM+^ (A 3 - v / 3A 8 )m+^m(a 3 - v^Ag) =0, (3.5a) 
1/(0, 0) = F(7r- 0,0), (3.5b) 
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(3.5c) 



for M = V and W as defined by Eqs. (3.3) and (3.4), with a trivial (0,0) 
dependence for V. 

Observe that the general matrix U (ip, fi, a, 6, 0) possesses the reflection sym- 
metries (3.5bc) only for -0 = and ir. The "point" ip — with U — V and its 
"antipode" ifj = n with U = W will be shown to correspond to the (unique) 
vacuum solution V and the new sphaleron S, respectively. As mentioned above, 
there is, in the radial gauge, a three-sphere with S at the top (tp = it) and V 
at the bottom (tp — 0), where the "height" corresponds to the energy of the 
configuration. 

The Cartesian components of the sphaleron- like SU (3) gauge field configu- 
ration are now defined as follows: 

A (r, 0, <t>) = , A m (r, 9, 0) = -g' 1 f(r) d m W(9, 0) W~ 1 (9, 0) , (3.6) 
with m = 1, 2, 3, and the following boundary conditions on the radial function: 

/(0) = 0, Km/(r) = l. (3.7) 
The corresponding Higgs field configuration is given by 



$(r,0,0) = h(r)r)W(9,(f)) 



I A 



v / 



(3.8) 
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with the following boundary conditions: 

h(0) = 0, lim/i(r) = l. (3.9) 

More specifically, the radial functions /(r) and /i(r) are even and odd in r, 
with /(r) oc r 2 and /i(r) oc r near the origin. 

For the vacuum solution V, the matrix W(9, 0) in Eqs. (3.6) and (3.8) is 
replaced by the constant matrix V defined by Eq. (3.3) and the profile function 
h(r) in Eq. (3.8) by the constant value 1. The vacuum field configurations are 
then A Q = A m = and $ = rj (1, 0, 0)\ with t standing for the transpose. 

The gauge and Higgs fields of Eqs. (3.6) and (3.8) are static and in the 
radial gauge (A r = 0). They carry the basic structure of the sphaleron S but 
are not general enough to solve the field equations consistently. In the next 
section, we present a suitable generalization. 



4 Sphaleron Ansatz 

4-1 Bosonic fields 

As shown in the previous section, the matrix function W(8, 0) from Eq. (3.4) 
has an axial symmetry given by Eq. (3.5a). This implies that the gauge field 
(3.6), with spherical components defined by 

A^-g- 1 f{r)d^WW~\ A e = -g' 1 f(r) d e W W' 1 , A r = , (4.1) 

has the property 



d^A c = [-2U 3 ,A C ], 



(4.2) 



where c stands for the component 0, 9, or r, and 2 U 3 is a constant matrix, 



-2U 3 = ~(X 3 - V3\ 8 ) =i 




1 
0-1 



(4.3) 



Next, introduce the following sets of matrices: 



T ( p = — sm — + cos <p — , 



Ai 



A 5 



T p = cos — + sin — , 



2i 



2i 



T3= 2i' 



(4.4a) 
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1/0 = + sin — + cos — , 



A, 



V p = cos — — sin — , 



2i 



2i 



3A 8 + A 3 



4i 



(4.4b) 



^ = sin(20) ^ + cos(20) ^ , C/ p = cos(20) ^ - sin(20) ^ , 



\/3 A 8 - A 3 
4i 



(4.4c) 



which have the same property 



d^X=[-2U 3 ,X], 



(4.5) 



with X standing for any of the matrices defined in Eqs. (4.4abc). The sets 
{T a }, {V a }, {U a } generate the three su(2) subalgebras of su(3) known as T- 
spin, V-spin, and [/-spin. Since T 3 , V3, and C/ 3 are not independent, we will use 
the basis {T^, T p , V^, V p , U^, U p , As/(2i), Ag/(2i)} of su(3) in the subsequent 
discussion. 

The gauge field (3.6) of the approximate sphaleron can be expanded in terms 
of this su(3) basis. It turns out that the component involves only T p , V p , 
Up, \s/(2i), Ag/(2i), and that Ag involves only T^, V^, U^. An Ansatz which 
has this structure and which generalizes (3.6) is given by 



gA (r,9,<f>) = 0, 



(4.6a) 



gA tj) (r,9,<j>) = a 1 (r,9) cos 9 T p + a 2 (r, 9) V p + a 3 (r,9) cos 9 U p + 

a A (r,9)^ + a 5 (r,9)^, (4.6b) 

g A e (r, 9, 0) = a 6 (r, 9) 7> + a 7 (r, 9) cos 9 + a$(r, 9) , (4.6c) 



gA r (r,9,<P) 



1 

r 



ag(r, 9) cos 9 + a w (r, 9) + a u (r, 9) cos 9 



, (4.6d) 



with real functions aij, for j = 1, . . . , 11, depending on r and 9 and matrices 
T, V, U from Eqs. (4.4abc) depending implicitly on 0. There are the following 
boundary conditions at the coordinate origin (r = 0): 



a j {0,6)=0, for j = 1,.. .,11, 



(4.7) 
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on the symmetry axis (9 = 0, it): 



aj(r, 9) = ctj(r) sin 9 
aj(r,9) = ctj{r) sin 2 # 
aj (r,9) = (-y~ 5 cos 9 d e a j_ 5 (r, 9) 
aj(r,9) = \ cos9de aj- 5 (r. 



for j = 1,2,9,10, (4.8a) 

for j = 3,4,5,11, (4.8b) 

for j = 6, 7, (4.8c) 

for j = 8, (4.8d) 



and towards infinity: 



lim 

r^oo 



c*i(r,0) ^ 




' -2 sin 0(1 + sin 2 #) 
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2 sin 9 cos 2 9 
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-2sin 2 # 


a 4 (r, 9) 




-sin 2 9(1 + 2 sin 2 fl) 
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sin 2 9 


a 6 (r,9) 
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017 (r, 9) 
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a 8 (r,9) 




-2 sin 9 


a 9 (r,9) 







a 10 (r,9) 







a u (r,9) j 








(4.9) 



Furthermore, the otj are required to have positive parity with respect to re- 
flection of the z-coordinate, 



aj(r, vr - 9) = +aj(r, 9), for j = 1, . . . , 11, 



(4.10) 



in order to keep the reflection symmetry of the approximate sphaleron config- 
uration; cf. Eqs. (3.5c) and (3.6). For completeness, the Cartesian components 



9 



of the gauge field Ansatz are 



9 A = - 



sin < 



P 

sin 9 cos <t> 



a\ cos9T p + a 2 V p + a 3 cos9U p + a 4 A 3 /(2i) + a 5 A 8 /(2i) 
cos #(a 6 + sin 9 a 9 ) + (cos 2 9 a 7 + sin 9 a 10 ) 



+ cos 9(a 8 + sin # an) 



(4.11a) 



2^2 



cos< 



a\ cos9T p + a 2 ^p + a 3 cos9U p + a4A 3 /(2i) + a 5 A 8 / (2i) 



+ 



sin 6* sin( 



cos 9 (chq + sin 9 a 9 ) + (cos 2 9 a 7 + sin # a 10 ) 



gA 3 = - 



+ cos 9(a s + sin a u ) U t 
cos 



sin a 6 — cos 9 a 9 ) + cos (sin a 7 — a 10 ) 



+ (sin 9 as — cos 2 an) Ua 



(4.11b) 



(4.11c) 



with a T-spin structure similar to the S and S* Ansatze [18]. 

Turning to the Higgs field (3.8) of the approximate sphaleron, we observe 
that the direction of the vacuum expectation value has been chosen so that 
2 ^^vac = 0. The field $ is found to have an axial symmetry, 



8^ = -2 u 3 $ , 



(4.12) 



which matches the axial symmetry (4.2) of the gauge field. The structure of 
this specific Higgs field is now generalized to the following Ansatz: 



$(r,0,0)=77 f3 1 (r,9)\ 3 + /3 2 (r,9) cos0 2iT p + (3 3 (r,9) 2i V p 



V / 



= 7] 



(3 2 (r,9) cos0e^ 
[P 3 (r,9)e- 1 * ) 



(4.13) 



with real functions (3k, for k = 1,2,3, depending on r and 9 and matrices 
T p and V p depending implicitly on <fi. The boundary conditions are at the 
coordinate origin (r = 0): 



de/3^0,9) = &(O,0) =/3 3 (O,0) = 0, 



(4.14) 
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on the symmetry axis {9 = 0, n): 



d e Pi{r,9) 



0, (3 k (r,9)=p k (r)sm9 



for k = 2,3, (4.15) 



and towards infinity: 



lim 



P2(r,9) 



( cos 2 ^ 

- sin# 

- sin# 



V 



(4.16) 



Furthermore, these functions must be even under reflection of the ^-coordinate, 

/3 k (r,Tr-9) = +/3 k (r,e), for k = 1,2, 3. (4.17) 

This completes the construction of the Ansatz for the bosonic fields. Be- 
fore discussing the resulting energy and field equations, the following three 
remarks may be helpful. First, the gauge field (3.6) and Higgs field (3.8) of 
the approximate sphaleron are reproduced by the following Ansatz functions: 







' -2 sin 0(1 + sin 2 #) ^ 


a 2 (r,9) 




2 sin 9 cos 2 9 


a 3 (r,9) 




-2 sin 2 9 


a 4 (r,9) 




- sin 2 9(1 + 2 sin 2 9) 


a 5 (r,9) 




VS sin 2 9 


a 6 (r,9) 


= f(r) 


2 


a 7 (r,9) 




2 
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a 10 (r,9) 
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(4.18a) 



fo(r,9) 
Ps(r,9) 



( cos 2 9 ^ 



h(r) 



— sin 9 

— sin 9 



(4.18b) 



which explains in part the choice of boundary conditions (4.7)-(4.9) and 
(4.14)-(4.16) for the general Ansatz. 
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Second, the Ansatz as given by Eqs. (4.6) and (4.13) has a residual 5*0(3) 
gauge symmetry under the following transformations: 



gA' n = Q (gA n + d n 



(4.19a) 



with 



exp 



u T (r, 9) 7> + u v (r, 9) V + uu(r, 9) 



(4.19b) 



for two-dimensional parameter functions u> s (r, 9), s — T,V, U. 

Third, the structure of the Ansatz gauge field (4.6) is quite intricate: for 
a given halfplane through the z-axis with azimuthal angle 0, the parallel 
components A r and A e involve only one particular su(2) subalgebra of su(3), 
whereas the orthogonal component A^ excites precisely the other five genera- 
tors of su(3). 



4-2 Energy and field equations 

The bosonic energy functional of theory (2.1) is given by 

E[A,Q]= f d 3 x -Itr(F mn ) 2 + Pm^r + A(|$| 2 -77 
Jr 3 v 



2 ^2 



(4.20) 



for spatial indices to, n = 1, 2, 3. From the sphaleron Ansdtze (4.6) and (4.13), 
one obtains 

A ^ /*oo /"^/^ 

E[A,^] = Air dr / d^ r 2 sin^ e(r,9) , (4.21) 

JO JO 

where the energy density e(r, 9) contains contributions from the Yang-Mills 
term, the kinetic Higgs term, and the Higgs potential in the energy functional, 

e(r, 9) = e YM {r, 9) + e Hk in(r, 9) + e Hpo t(r, 9) . (4.22) 

The detailed expressions for these energy density contributions are relegated 
to Appendix A. The energy density (4.22) turns out to be well-behaved due 
to the boundary conditions (4.7)-(4.9) and (4.14)-(4.16). In addition, there 
is a reflection symmetry, e(r, 9) = e(r, n — 9), which allows the range of 9 in 
Eq. (4.21) to be restricted to [0,tt/2]. 

It can be verified that the Ansdtze (4.6)-(4.10) and (4.13)-(4.17) consis- 
tently solve the field equations of the theory (2). That is, they reproduce the 
variational equations from the Ansatz energy functional (4.21). This result 
is a manifestation of the principle of symmetric criticality [25], which states 
that in the quest of stationary points it suffices, under certain conditions, to 
consider variations that respect the symmetries of the Ansatz (here, rotation 
and reflection symmetries). 
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It is difficult to prove that the resulting system of partial differential equa- 
tions for the Ansatz functions Oij(r,9) and /3k(r,0), together with the condi- 
tions (4.7)-(4.10) and (4.14)-(4.17), has a nontrivial solution, that is, a non- 
vacuum solution. It is, however, possible to give a heuristic argument in favor 
of a nontrivial regular solution (a similar argument has been used for S* [18]). 

According to the boundary condition (4.16), the function f3\(r,9) vanishes 
asymptotically for 9 = n/2, i.e., on the equatorial circle at infinity. Let us 
simply exclude the case of an "isolated" zero of f3i(r,6) at (r,9) = (oo,7r/2), 
which is discussed further in Appendix B. For a regular solution and by reflec- 
tion symmetry 6 — > n — 9, the set of zeroes of (3i(r, 9) then extends inwards 
from infinity along the equatorial plane 9 = tt/2 and intersects the symmetry 
axis of the Ansatz at the coordinate origin r = 0. Since /^(r, 9) and (3 3 (r,9) 
also vanish due to the boundary conditions (4.15), the Higgs field (4.13) is 
exactly zero at this point, $(0,0,0) = 0, which is not possible for vacuum 
configurations with |<&(r, 9,(f>)\ = rj. 

It appears that the details of the (regular) solution can only be determined 
by a numerical evaluation of the reduced field equations. The explicit numer- 
ical solution of these partial differential equations is, however, rather difficult. 
In the following, we will simply use the approximate sphaleron fields (3.6) and 
(3.8), which are sufficient for our purpose of looking for fermion zero modes. 2 



5 Fermion zero modes 

5.1 Ansatz and zeromode equation 

As mentioned in Section 2, the SU(3) gauge theory considered has a single 
massless left-handed fermionic field in the 3 representation. This fermionic 
field is now taken to interact with the fixed background gauge field of the 
approximate sphaleron (3.6) or the general sphaleron Ansatz (4.6). In both 
cases, the corresponding fermion Hamiltonian, 

H = ia-(V-gA), (5.1) 

has the following symmetries: 

[K 3 ,H}=0, {R 3 ,H} = 0, (5.2) 



2 The simplified fields (3.6) and (3.8), with appropriate radial functions f(r) and 
h(r), can also be used to get an upper bound on the energy of S. For the case 
of X/g 2 = 0, one obtains E[S] < 1.72 xB[S], where E[S] denotes the energy of 
the SU(2) sphaleron S [11] embedded in the SU(3) Yang-Mills-Higgs theory. This 
implies that E[S] < E[S*], at least for X/g 2 = 0, where E[S*] « 1.91 x E[S] 
denotes the energy of the embedded SU(2) sphaleron S* [18]. 
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with 



>^3 , 



R 3 = diag(-l, +1, -1, +1, -1, +1) TZ 3 , 



(5.3a) 
(5.3b) 



and IZ3 the coordinate reflection operator with respect to the 3-axis, i.e., 
IZ3 9 = 7v — 9, TZ 3 4> = 0, and TZ 3 r = r. The A and a matrices in Eq. (5.3a) 
operate, of course, on different indices of the fermionic field, the SU (3) "color" 
and SU{2) spin indices respectively (see below). Observe that R3 and K 3 
commute, [R3, K 3 ] = 0. 

The eigenvalue k of K 3 can take odd-half-integer values (k e Z + \) and 
the corresponding eigenstate has the form 



§M,0) 



F R+ (r,9 
F G+ (r,9 
F B+ (r,9 
F R -(r,9 
F G -(r,9 
F B -(r,9 



e i(«-l/2)4. 
e i(/c+l/2)<^ 
gi(«-3/2)4> 
gi(K+l/2)</> 
e i (k+3/2) 
p i(«-l/2)4> 



(5.4) 



where + and — denote the spin (in units of h/2) and R, G, and B the color 
(red, green, and blue). The real functions F cs (r, 9) must vanish for r — * 00, in 
order to have a normalizable solution, 



dr d9 d<p r 2 sin9 |^| 2 = 1. 

JO JO JO 



(5.5) 



Those functions F cs (r, 9) which are accompanied by a ^-dependent phase fac- 
tor in Eq. (5.4) are required to vanish also for 9 = 0, 7r and r = 0, in order to 
have regular field configurations. 

For the gauge field (3.6) of the approximate sphaleron and a particular 
eigenvalue k of the Ansatz (5.4), the zeromode equation H^/ = reduces to 
the following partial differential equation: 

{MO) ^ + \ M9) ^ + \ A^) + j -^L Ai{e) I *( r , 0, o) = , (5.6) 

with matrices 



MO) = 



( 



sin6> cos 9 I3 
sin 2 9 1 3 



sin 2 1, 



sin6> cos 9 1 3 



(5.7a) 
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A 2 (9) = 



' -sin 2 01 3 sin# cos#l 3 , N 
^sin#cos#l 3 sin 2 9 1 3 



(5.7b) 



A 3 (n) = 







-re + 1/2 










-re- 1/2 












re + 1/2 
re + 3/2 









-re + 3/2 













re - 1/2 






, (5.7c) 



A A {9) 



( aW 



AT (9) 



A?\e) 



-Af\6) A { 1 ] \6) ) 



(5.7d) 



4 1} (0) 



sin 2 # cos 9 sin 2 9^ 

- sin 2 9 sin 3 9 

- sin 2 9 cos 9 sin 3 6> 



(5.7e) 



A?\9) 



sin 4 # 



sin 3 9 cos # —2 sin 9 cos 2 # \ 



sin 9 cos 0(3- cos 2 9) sin 2 (cos 2 0-2) 2 cos sin 2 9 



V 



sin 2 



(5.7f) 



where the superscript t in (5.7d) indicates the transpose. Fermion zero modes 
correspond to normalizable solutions of Eq. (5.6). 

Since R 3 anticommutes with H, {i? 3 , H } = 0, common eigenstates with 
zero energy can be found. As shown in Appendix C, zero-energy eigenstates 
with opposite eigenvalues of K 3 and i? 3 can be constructed if one starts from 
a generator U of n 5 [SU(3)] with opposite winding number compared to U. 

5.2 Asymptotic behavior 

The gauge field function /(r) of the sphaleron-like field (3.6) is approxi- 
mately equal to one for large enough r. With / = 1 exactly, the Yang-Mills 



15 



field is pure gauge, 



A m = -g- 1 d m WW- 1 , 



(5.8) 



in terms of the SU{3) matrix W(9, <p) defined by Eq. (3.4). The Dirac Hamil- 
tonian becomes then 



H = ia-(y-gA) = WH ircc W~\ 



(5.9) 



with H{ rcc = i<7 • V. Asymptotically, the normalizable zero mode can therefore 
be written as 



&°°\r,e,<f>) = w(e,cf>) * free (r, 0, 0) , 



(5.10) 



for an appropriate solution of H{ rce \l/f r ee = 0. 

The wave function $( 00 ) is an eigenstate of K 3 , if \I>f ree is an eigenstate of 
K 3 with the same eigenvalue, where K 3 is defined by K 3 = — id^ — \{\$ — 
As) + \(Tz- Since H{ rcc does not mix colors, the asymptotic two-spinors can 
be given separately for each color, 



41(^,0) 



„i+i 



+ v /(/ + l/2- K -c)/(2/ + l) ^ +c - 1/2 (^,0)\ 



^ + 1/2 + k + c)/(2l + 1) lf +c+1 / 2 (0, 0) ; 



. (5.11) 



Here, k is the ^-eigenvalue, / stands for the orbital angular-momentum quan- 
tum number, and c takes the values 0, —1, +1 for the colors R, G, B, respec- 
tively. Most importantly, the quantum number I is restricted by the condition 



-1<K + C-±<1-1, 



(5.12) 



which ensures the consistency of the orbital angular momentum and the 0- 
dependence from the i^ 3 -symmetry. 

The asymptotic i^ 3 -eigenstate ^(°°)(r, 9, <fi) with eigenvalue k is then given 
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by 



W 



\ 



E 
E 
E 
E 
E 
E 



oc 



oo 



k| + t 



D (l) f 

£>g) r -(J+l) / 

D m r - { i +1) f 

lie+ll+5 B V 



Z + ±-/c)/(2Z + l) If- 1 /^ 

;z + |-«)/(2z + i) i?- 3/2 



/ - | - k)/(2Z + 1) >f +1/2 



Z + i + «)/(2Z + l) lf +1/2 



Z-| + «)/(2Z + 1) l?" 1 ' 2 



Z + | + k)/(2Z + 1) Y ; K+3/2 



(5.13) 



with constant coefficients , Dq 1 , and -D^. This implies an asymptotic r~ 2 
behavior for the fermion zero modes, provided the sums in (5.13) converge 
and are nonvanishing. 

For r close to zero, one has f(r) oc r 2 from the requirement of finite Yang- 
Mills energy density. With / = exactly, the gauge field vanishes, A m = 0. 
The behavior at the origin is then given by 



$(°) 



v^oo 

^1= 



^1= 



E 
E 
E 



V 







/(Z + I + «)/(2Z + 1) ^" 1/2 + 4 0) <^-±,o 




^\ 


/(Z + | + «)/(2Z + l) ^ +1/2 + 4 0) ^,o 






/(Z - \ + k)/(21 + 1) YT 3/2 + 4 0) ^-|,o 


i+l 


^\ 


/(Z + I - «)/(2Z + 1) ^ K+1/2 + 4 0) ^ + i,o 




^\ 


/(/ - I - «)/(2Z + 1) ^ +3/2 + 4 0) <W§ ,o 






/(Z + § - «)/(2i + 1) YT 1/2 + 4 0) <^-±,o 



(5.14) 

with constant coefficients 4\ 4^ an d 4^- The wave function \J/(°)(r, 0, 0) 
is a i^3-eigenstate with eigenvalue k. 

An explicit solution of the zeromode equation (5.6) would relate the coef- 
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ficients D and C of Eqs. (5.13) and (5.14). An analytic solution is, however, 
not feasible, possibly except on the symmetry axis. 



5.3 Solutions on the symmetry axis 

The zeromode equations can be solved on the symmetry axis (coordinate 
z = r cos 9) for the background gauge field (3.6) of the approximate sphaleron. 
Two zero modes have been found, one with k — +1/2 and one with k = —1/2. 

For k = — 1/2, we can, in fact, construct an analytic solution. On the z-axis, 
only the ^-independent components of the general Ansatz (5.4) are nonzero. 
This leads to the simplified Ansatz 



# (_) (0,0,. 



1 

72 



lt\z) 








(5.15) 



for cylindrical coordinates (p, 0, z) and with even functions 4 \ z ) and l{ \z). 
The wave function is an _R 3 -eigenstate with eigenvalue 1. Inserting (5.15) 
into the zeromode equation (5.6) with all 9 derivatives set to zero, one obtains 
the normalizable solution 



l(z), 



(5.16) 



in terms of the following function of z: 



iW _ P (-/>^>). 



(5.17) 



Hence, |$(~^| is normalized to one at the origin and drops off like \z\~ 2 for 
large \z\. Note that the same exponential factor appears in the fermion zero 
mode of the spherically symmetric SU (2) sphaleron S [13-17]. Similar analytic 
solutions of fermion zero modes on the symmetry axis have also been found 
for the axially symmetric SU(2) constrained instanton I* [19,20]. 
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Fig. 3. Left panel: radial function f(Q from Eq. (5.19) for the SU(3) back- 
ground gauge field (3.6), with /(C) shown as a function of the compact coordinate 
£ = r/(r + 1). Right panel: functions of the corresponding k = —1/2 fermion zero 



mode (5.15) on the symmetry axis (r 
broken curve. 



z\), with 4 ^(C) = ^4 ; (C) shown as the 



For k — +1/2 and 9 close to or ir, we use the Ansatz 

(l{ + \r) cos 6 \ 
4 +) (r) sinfle^ 
4 +) (r) sin# cos0e _i * 



$(+)(r,0,0) = 



4 +) (r) sinfle^ 
4 +) (r) sin 2 cos#e 2i ^ 



(5.18) 



The wave function is an i? 3 -eigenstate with eigenvalue 1 which is nor- 
malized to one at the origin, provided 4 +) (0) = 1 and Z«(0) = for 
n — 1, . . . , 5. The differential equation (5.6) can then be solved numerically 
to first and second order in 6 or 7r — 6. On the symmetry axis, only the first 



and last entries of (5.18) are nonzero, but the functions 4 + \ • • • , s ^ a f" 



;(+) 



feet the equations for Z^ and l^' because of the coupling through the 8 
derivatives in (5.6). On the z-axis, one has asymptotically oc \z\~ 2 for 

(Zi,Z 2 ,Z3,Z4,Z 5 ,Z 6 )cx(-l,l,0,0,0,l)/r 2 . 
In order to be specific, we take the following radial gauge field function: 



/(+) 



/(C) 



(4C-l)7[(4C-l) 2 + (4C-3) 2 ] 



if C < 1/4 

if Ce (1/4,3/4) , (5.19) 
if C > 3/4 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Fig. 4. Numerical solutions for the functions of the k = +1/2 fermion zero mode 
(5.18) on the symmetry axis, with SU(3) background gauge field given by Eqs. (3.6) 
and (5.19) for £ = r/(r+ 1) and r = \z\. Left panel: 4 + ^(C) an d ~~ 4 + ^(0 are shown 
as broken and solid curves, respectively. Right panel: 4 + ^(C) an d 4 + ^(0 are shown 
as broken and solid curves, respectively [not shown are 4 + '*(C) an d 4 + ^(C)i which 
take values close to zero]. 

in terms of the compact coordinate ( G [0, 1] denned by 

C = r/(r+l). (5.20) 

The numerical solutions for the functions of the k = —1/2 and k = +1/2 
zero modes on the axis are shown in Figs. 3 and 4, respectively. A complete 
numerical solution of these zero modes over the (r, #)-halfplane is left to a 
future publication. 

The analysis of this section does not rigorously prove the existence of fer- 
mion zero modes. It is, in principle, possible that the only solution of the 
zeromode equation (5.6) which matches the asymptotic behavior (5.13) is the 
trivial solution vanishing everywhere. But the existence of two normalizable 
solutions on the symmetry axis, which is already a nontrivial result, suggests 
that the partial differential equation (5.6) has indeed two normalizable so- 
lutions. In fact, the symmetries of the S Ansatz may play a crucial role for 
the existence of fermion zero modes, just as was the case for the S and S* 
sphalerons. 



6 Conclusion 

As explained in the Introduction and in more detail in Ref. [8], one pair of 
fermion zero modes (with a cone-like pattern of level crossings) is needed to 
recover the SU(3) Bardeen anomaly as a Berry phase (Fig. 2). The existence 
of one pair of fermion zero modes somewhere inside the basic noncontractible 
sphere of gauge-transformed vacua follows from the family index theorem. In 
this article, we have argued that these zero modes are, most likely, related to 
a new sphaleron solution of SU(3) Yang-Mills-Higgs theory, indicated by S 
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in Fig. I. 

The self-consistent Ansdtze for the bosonic and fermionic fields of S have 
been presented in Sees. 4.1 and 5.1. Specifically, the gauge and Higgs fields 
are given by Eqs. (4.6)-(4.10) and (4.13)-(4.17), and the fermionic field by 
Eq. (5.4). It appears that the reduced field equations can only be solved nu- 
merically. The sphaleron S has at least three negative modes, if it is indeed 
at the "top" of a noncontractible three-sphere in configuration space [cf. the 
discussion in the paragraph starting a few lines below Eq. (3.5)]. 

In pure SU (3) Yang-Mills theory, there exists perhaps a corresponding Eu- 
clidean solution, just as the BPST instanton I [26] corresponds to the sphaleron 
S (the gauge field of a three-dimensional slice through the center of I roughly 
matches that of S). This non-self-dual solution would have at least two neg- 
ative modes. It is, however, not guaranteed that the reduced field equations 
have a localized solution in IR 4 . For a brief discussion of non-self-dual solutions, 
see, e.g., Section 6 of Ref. [27] and references therein. 

Purely theoretically, it is of interest to have found the sphaleron S in SU (3) 
Yang-Mills-Higgs theory which may be "responsible" for the non-Abelian 
(Bardeen) anomaly. More phenomenologically, there is the possibility that S- 
like configurations take part in the dynamics of Quantum Chromodynamics 
(QCD) or even grand-unified theories. The SU (3) Bardeen anomaly cancels, 
of course, between the left- and right-handed quarks of QCD, but the config- 
uration space of the gluon gauge field still has nontrivial topology. Quantum 
effects may then balance the SU(3) gauge field configuration given by Eqs. 
(4.6)-(4.10). The resulting effective sphaleron (Sqcd) can be expected to play 
a role in the nonperturbative dynamics of the strong interactions. 



A Ansatz energy density 

In this appendix, the energy density from the Ansatz gauge and Higgs fields 
of Eqs. (4.6) and (4.13) is given. Specifically, the contributions to the energy 
density (4.22) from the Yang-Mills term, the kinetic Higgs term, and the Higgs 
potential term in the energy (4.20) are 



cym — 



2 g 2 r 2 sin 2 9 



cos 2 9 d r cti + (ctg + O4O9 — 7jOt20tn — ^03010) h 
+ (9 r «2 — (oiio — lectio + I cos 2 9 0:30:9 — r\V3a 5 a w — \ cos 2 9 a^nj /? 

<9r03 — (20n — |o 2 Og + |o 4 On — Iv^OsOii — TjOiOio) /r 

d r a>4 + (| cos 2 9 03O11 — cos 2 9 01O9 — |o 2 oio) / r 
+ d r a 5 — ^VS (o 2 o 10 + cos 2 9 a 3 a u j /■ 



+ cos 2 9 
+ 
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+ 

+ 
+ 

+ 

+ 
+ 
+ 
+ 



2g 2 r 2 



d r a> 6 — (cos 9 dgag — \olsOl w + \ cos 2 9 a 7 « n — sin 9 o 9 ) / r 



cos 9 d r a 7 — (dgaio — \ cos 9 o 6 On + | cos 9 o 8 o 9 ) jr 
d r a>g — (cos 9 doa>n + ia 6 oio — \ cos 2 9 a-jag — sin 9 a±ij jr 



1 



2 g 2 r 4 sin 2 9 



«6 — sin 9 «i — 4o 2 o 8 + O4O6 — \ cos 2 6? 030:7 + cos 9 dgai 



l2 



cos 9 o 7 + i cos # (03O6 — oio 8 — a/3 0:507 — 04O7) — <%o 2 
2o 8 + sin #03 + |o 4 o 8 — |o 2 o 6 — |v^3 05O8 — \ cos 2 # o x o 7 — cos 9 dea 3 
cos 9 (aids + \a 2 a 7 — |o 3 o 8 ) — dga 4 

cos9 (o 3 o 8 + o 2 o 7 ) — ^05 1, (A.l) 



eHkin = V 



d r Pi — \ [cos 2 9a 9 f3 2 + l r + cos 2 # 



X 



drfo + \ (ag/3i ~ aii/3 3 ) A + <9 r/ 3 3 + \ ( «io/?i + cos 2 9 «nft) /r 



— 



dePi - I cos 9 (o 7 /3 3 + o 6/ 3 2 ) 1 + [ d e j3 3 + \ cos 9 (o 8 /3 2 + o 7/ 9i) 



n 2 



+ 
+ 



cos 9 d e (3 2 - sin 9 f3 2 + \ (o 6/ 9i - o 8/ 9 3 ) 



77^ 



r 2 sin 2 9 1 4 



o 4( #i + (l/\/3 ) o 5 /3i + cos 2 9 ax(3 2 + o 2( 5 3 



+ \ cos 2 # 2/3 2 - oi/3i + o 4 /3 2 - (l/\/3 ) o 5/ 3 2 - o 3( 5 3 
+ i [ 2(3 3 + o 2 A - (2/ y/3 ) o 5/ 9 3 + cos 2 9 a 3 f3 2 ] 2 J , 



e Hp ot= Ar? 4 /3 2 + C os^^ + /3|-l 



2 /) <?2 



(A.2) 



(A.3) 



These energy densities are manifestly symmetric under reflection 9 — > 7r — 
provided the Unsafe functions Oj(r, #) and j3k(r,9) are; cf. Eqs. (4.10) and 
(4-17). 



B Argument for a nontrivial regular solution 

In this appendix, we present a heuristic argument for the existence of a 
nontrivial regular solution of the reduced bosonic field equations, with Ansatz 
functions aj(r,9) and (3 k (r,9) from Eqs. (4.6)-(4.10) and (4.13)-(4.17). The 
reasoning here complements the one of Section 4.2. It focuses on different 



22 



energy contributions from the region near the equatorial plane [0 = tt/2 or 
z = 0] and uses a simple scaling argument to show that the total energy of 
a regular solution cannot be zero. The coupling constant \/g 2 is taken to be 
nonzero. 

First, consider the possibility that (3i(r,9) has a vacuum value at the co- 
ordinate origin, /3i(0,7r/2) = ± 1. We now try to get a vacuum solution over 
the whole (r, #)-halfplane but do not allow for singular configurations with a 
finite change of /3i(0, tt/2) over an infinitesimal r interval at the origin r = 0. A 
vanishing Higgs potential energy density (A. 3) on the equatorial plane then re- 
quires /3 3 (r, 7r/2) = — [ 1 — Pi(r, tt/2) 2 } 1 / 2 , assuming (3\(r, tt/2) to change mono- 
tonically from the value ± 1 at r = to the asymptotic value as given by the 
boundary condition (4.16). Define the length scale R by [3s(R, tt/2) = — 1/2 . 

With 01 + /3| ~ 1 near the equatorial plane, a vanishing first curly bracket 
in the kinetic Higgs energy density (A. 2) requires a w ~ —2(rd r (3 3 )/fli ~ 
2(rd r (3i)/(3 3 , which peaks somewhere around r = R. The corresponding 
contribution of a±o to the Yang-Mills energy typically scales as \/R. This 
energy contribution can be reduced by increasing the value of R. But for 
large R and /3 k (r,9) of the asymptotic form (4.18b), the ^-dependence of 
—2(rd r p 3 )/Pi and 2(rd r (3i)/fl 3 is entirely different. Hence, the previous ex- 
pression for a lo (r,0) needs to be modified near 9 tt/2 for large r and the 
kinetic Higgs energy term picks up a contribution which typically scales as R. 
Clearly, both terms cannot be reduced to zero simultaneously by changing R. 
This basically rules out having a regular vacuum solution with |/3i(0, tt/2)\ = 1. 

Second, consider the possibility that the value of \Pi(0, tt/2) \ does not equal 
1. Since fair, 9) and /5 3 (r, 6) vanish on the whole z-axis by the boundary condi- 
tions (4.15), one then has for the Higgs field (4.13) at the origin |$(0, 0, 0)| ^ rj, 
which already differs from the vacuum solution with |$(r, 9,(f>)\ = r\ every- 
where. 

Apparently, both possibilities lead to a regular solution with nonzero total 
energy, which concludes our heuristic argument. The actual solution may very 
well have a Higgs field $(r, 9, (ft) that vanishes at the coordinate origin, with 
a corresponding localized nonzero energy density; cf. Section 4.2. 



C Fermion zero modes from gauge fields with winding number — 1 

Instead of the generator U of tt 5 [SU(3)] with winding number +1 [given by 
Eq. (3.1)], one may use a generator U with winding number —1. This matrix 
U is simply the inverse of U, namely U(ip, /i,a,0,(f)) = f/~ 1 (-^, fi,a,9,(p). The 
background gauge field for the fermion zero modes is then determined by the 
matrix function W(9,(p) = W~ 1 (9,(f)), with gA^ = — f(r) d^W W' 1 for the 
approximate anti-sphaleron. 

The corresponding Dirac Hamiltonian H = i a m (d m — gA m ) has the sym- 
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metry properties 

[K 3 ,H}=0, {R 3 ,H} = 0, 



(C.l) 



(C.2) 



where 

K 3 = -idf - \ (A 3 - A 8 ) + |cr 3 , 

and i?3 is still given by (5.3b). An eigenstate of K 3 with eigenvalue k has the 
form 

/ tt* /_ /IN _i Ck-I /21<A \ 



F R+ (r,e) 


e i(K-l/2)«*> 


F G+ (r,9) 


e i (n-3/2) 


F B+ (r,0) 


e i(K+l/2)<^ 


F R -(r,0) 


e i(K+l/2)<^ 


F G -(r,6) 


e i («-l/2) 


F B -(r,6) 


gi (k+3/2) 



(C.3) 



The zeromode equation H iff — now reduces to 

^A 1 (e)-^ + ^A 2 (9)^ + ^A 3 (K) +Ma 4 (0)} *(r,fl,0) = 0. (C.4) 

Here, the matrices Ai and A 2 are given by Eqs. (5.7a) and (5.7b), the matrices 
A 3 and A 4 by 

A 3 {k) = NA 3 {-k) AT 1 , I 4 (0) =A^A 4 (7r-^)A^- 1 , (C.5) 
with A 3 and A 4 from Eqs. (5.7c) and (5.7d), and the constant matrix iV by 



N 



I 








(i 





0^ 








o 3 







-1 









fl 





o\ 


v° 





- l ) 









-1 







o 3 






V 


1° 





-v 








) 



(C.6) 



Note that the matrices A 1 and A 2 have the following reflection properties: 
A 1 (6) = N A x {ti - 6) N~\ A 2 {9) = -N A 2 (n - 6) N- 1 . (C.7) 
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Hence, there is a one-to-one correspondence between the solutions of Eqs. (5.6) 
and (C.4): \&(r, 0,0) solves (5.6) and has -/^-eigenvalue k if and only if 
\I/(r, 9, 0) = N^(r, n — 9, 0) solves (C.4) and has ^-eigenvalue —k. 
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